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ABSTRACT 

The vast diversity of planetary systems detected to date is defying our capability of understanding their formation and evolution. 
Well-defined volume-limited surveys are the best tool at our disposal to tackle the problem, via the acquisition of robust statistics of 
the orbital elements. We are using the HARPS spectrograph to conduct our survey of as 850 nearby solar-type stars, and in the course 
of the past nine years we have monitored the radial velocity of HD103774, HD109271, and BD-061339. In this work we present 
the detection of five planets orbiting these stars, with m sin(/) between 0.6 and 7 Neptune masses, four of which are in two multiple 
systems, comprising one super-Earth and one planet within the habitable zone of a late-type dwarf. Although for strategic reasons 
we chose efficiency over precision in this survey, we have the capability to detect planets down to the Neptune and super-Earth mass 
range, as well as multiple systems, provided that enough data points are made available. 

Key words, extra-solar planets - multiple planetary systems - radial velocity - HD103774 - HD109271 - BD-061339 



1. Introduction 



After more than 15 years of exo-planets detections ( Mayor 1995 
ISchneider 20121 1. we have witnessed an extreme diversity of 
planetary systems, but we have not yet detected a system that 
could be considered a twin to our own solar system, not even 
limited to the giant (and easier to detect) planets. We are left 
with the question of whether our system, and with it our planet, 
are the norm, or rather the exception in planet formation and evo- 
lution. To date, only about one third of the detected planets are 
found in confirmed multi-planet systems (Schnei der 2012b . This 
fraction is steadily rising with the increasing amount of available 
high-quality data, indicating that we are fighting against a detec- 
tion bias. Multi-planet systems are more challenging to detect 
because they require many data points to permit the identifica- 
tion of the individual signals and because, as is revealed from 
the data collected so far, often they involve low-mass planets 
in the super-Earth to Neptune mass range (e.g. ILo Curto 2010 
IHebrard 20101 ILovis 20TT1 IFischer 20121 >. The combination of 
the low amplitude of these signals with the more complex mo- 
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tion induced on the parent star by the systems requires a substan- 
tial observational effort both in terms of observing time and of 
the measurement precision. 

Our survey is using the HARPS spectrograph, which is 
mounted full time at the 3.6m ESO telescope in La Silla 
and is optimized for precise radial velocity measurements 
( |Mayor 2003| >. To date, it has discovered more than 150 exo- 
planets, most of them (through the high-precision program) with 
masses in the super-Earth range; this definitely changed our view 
on the field. 

Since the start of HARPS operations we are monitoring 
about 850 stars with moderate radial velocity precision (« 2 — 
3ms~ l ), searching for Jupiter and Neptune planets orbiting solar- 
type stars (Lo Curto 2010), while the HARPS high-precision 
program (Mayor 2003), which routinely achieves precision bet- 
ter than lms" 1 , targets lower mass planets, mostly in the super- 
Earth mass regime. The goals of our program are the acquisi- 
tion of accurate orbital elements of Jupiter and Neptune-mass 
planets in the solar neighborhood (within 57.5pc) and the search 
for Jupiter twins, which might reveal solar system twins. Of the 
43 planets discovered so far by this program, 17 have periods 
longer than 1000 days dMoutou 20091 INaef 20101 two more are 
presented in a paper in preparation), but their eccentricities are 
significantly higher than that of Jupiter. The program has de- 
tected three systems with two planets and one triple planetary 
system so far (IHebrard 20101 ILo Curto 20101) . Out of the nine 
planets in multiple systems we have discovered so far, only two 
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have masses larger than or equal to Jupiter, which is consistent 
with the observed trend of low-mass planets being detected pri- 
marily in multiple systems. 

In this paper we report three stars with multiple Keplerian- 
like signals. In section 2 we discuss the stars' properties, and in 
section 3 the orbital solutions. In section 4 we briefly discuss the 
encountered systems. 

2. The stars 

The three stars belong to our volume-limited sample and 
were selected to be low-activity, non-rotating, solar-type stars 
on the main sequence dLo Curto 201 01) . We have used the 
magnitude estimates by Hog (2000), and the one by Koen 
(2011) for the redder star BD-061339. They originate from 
the Tycho and HIPP ARCOS catalogs, respectively dESA- 19971 
IVan L eeuwen 2007). The color indices and the parallaxes are 
from lVan Leeuwen 20071 and the bolometric corrections are ex- 
tracted from a calibration to the data presented by Flower (1996), 
with uncertainties generally below 0.1 mag. 

2.1. HD103774 

The star HD 103774 is a young, metal-rich, main-sequence star 
with spectral type F5V (Pic kles 20101 ISousa 201 II) . Its average 
activity level is low, but it varied considerably during the time 
span of our observations (see figure 0), going from a more "ac- 
tive" state with logiRnK) — -4.72 to a more "quiet" state with 
log(R HK ) = -5.00 (see |Noyes 1984| for a definition of the ac- 
tivity indicator). We therefore expect to measure some level of 
radial velocity jitter. 

Using the online data service VizieR, we found measure- 
ments of this star in the mid- and far-infrared, from the IRC (in- 
strument onbo ard the AKARI satellite, Ishih ara 20T0l l and IRAS 
dMoshir 19 90) catalogs (see figure [T}. The data points at 25, 60, 
and 100 pmx are upper limits. The solid line in figure[T]reproduces 
the black body corresponding to the star's photospheric tempera- 
ture of 6489/T. The uncertainty on the temperature (77K) would 
not be noticed on the scale of this graph. The point at 12/j.m 
is deviating from the black body and might seem to indicate 
an infrared excess. Circum-stellar dust is indeed quite common 
around early-type young stars. While the evidence for infrared 
excess at 12 /jm is only marginal (1 sigma) and calls for more 
precise measurements, no excess is measured up to 9 fim, giv- 
ing an indication on the minimum size of the grains that could 
possibly be present in the disk. 

2.2. HD109271 

This old and quiet G5V star does not show infrared excess 
(Ishihara 2010 IRAS data are not available), and has a near- 
solar metallicity (ISousa 201 II) . Its basic parameters, as for the 
other stars, are shown in table [T] 

2.3. BD-061339 

BD-061339 is a late-type dw arf (spectral ty pe K7V/M0V), with 
no obvious infrared excess (Mos hir 19901 IRAS Faint Source 
Reject Catalog, Version 2.0). 

The estimation of the metallicity in late-type dwarfs is com- 
plicated by the abundance of molecular bands of various metal 
oxides in the spectra (e.g. TiO, VO). Neves et al. (2012) have 
discussed various calibrations present in the literature and pre- 
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Fig. 1. Spectral flux density for HD 103774 measured in Jy. The 
solid line is a black body with the temperature of the star: 6489K. 
The arrows indicate that the three data points are only upper lim- 
its. 

sented an improved photometric calibration for the metallicity. 
We adopted their prescription to estimate the metallicity of BD- 
061339, and used the dispersion around their calibration as an 
estimate of the uncertainty. This value is well above the mea- 
surement uncertainty computed via error propagation from their 
formula, and is limited by the quality of the calibration. 

Owing to the red color of the star the estimate of the activity 
index log R' UK is considered unreliable and no value is included 
in table Q] 

3. Orbital solutions 

We studied the radial velocity (RV) variation of our target stars 
via the analysis of the periodograms of the RV data following 
the generalized Lomb-Scargle periodogram method outlined by 
Zechmeister and Kuerster (2009) and via the use of genetic al- 
gorithm as described in the work by Segransan et al. (201 1). In 
all cases we verified whether the RV was correlated with the bi- 
sector span of the cross correlation function (CCF), the activity 
indicator R HK ( |Noyes 1984) , and the FHWM of the CCF be- 
fore claiming a companion. In addition to this, we also studied 
the periodogram of these indicators. When attempting a fit with 
more degrees of freedom (e.g., using a model with more planets 
or adding a long-term drift to the model), we verified with an 
F-test that the new fit brings a significant improvement over the 
fit with fewer degrees of freedom. We ran an F-test and obtained 
a false-alarm probability (FAP) for which we used a threshold 
of 10~ 4 to identify significant peaks. The orbital elements are 
finally displayed in table [2] 

3.1. HD1 03774 

We acquired 103 radial velocity measurements for this star in 
about 7.5 years. The average uncertainty is ss 2.6ms~ l , while the 
standard deviation of the radial velocity over the entire data set is 
* 26m.iT 1 , a clear indication of variability. The periodogram of 
the radial velocities (Pig. |2j indicates an excess of power at » 5.9 
days, while no other significant peak is visible. The one-year and 
the two-year aliases (at 5.79d and 5.84d) are also present in the 
periodogram, but at half the power of the main signal. 

In figure [3] we show the Keplerian fit as a function of phase 
together with the data. The orbit is evident, although noisy. 

Although the periodogram shows various structures 
("bumps" at * lid, 13d, 15d, 470d, 740d), they are not signifi- 
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Table 1. Observed and inferred stellar parameters of the three stars of this study, n: parallax, d: distance, My: absolute magnitude in 
V, Bol. corr.: bolometric correction, L: luminosity, T e ff. effective stellar temperature, M*: stellar mass, log g : gravity at the star's 
surface, vsin i: projected rotational velocity of the star, log R' RK : activity indicator fNoyes 1984| l, f rot (log ^ K ): star's rotation period 
estimated via the activity indicator. 



Parameter 




HD 103774 


HD 109271 


BD-061339 




Hipparcos name 


HIP 


58263 


61300 


27803 




Spectral type 




F5V 


G5V 


K7V / MOV 




V 


[mag] 


7.12 ± 0.01 


8.05 ± 0.01 


9.69 ± 0.01 


dHog 2000||Koen 2010b 


B- V 


[mag] 


0.503 ± 0.014 


0.658±0.002 


1.321 ±0.001 


Jvan Leeuwen 20071 


n 


[mas] 


18.03±0.52 


16.21±0.75 


49.23±1.65 


iVan Leeuwen 2007b 


d 


[pc] 


55±2 


62±3 


20±1 




My 


[mag] 


3.40 ± 0.06 


4.1 ± 0.1 


8.15 ±0.07 




Bol. Corr. 


[mag] 


-0.013 


-0.097 


-0.782 


JFlower 1996b 


L 


Ho] 


3.5 ±0.3 


2.0 ±0.3 


0.095 ± 0.01 




T e g 


[K] 


6489 ± 77 


5783 ± 62 


4324 ± 100 


JSousa 201 lllAmmons 2006b 


[Fe/H] 


[dex] 


0.28 ± 0.06 


0.10 ±0.05 


-0.14±0.17 


(Sousa 201 1 , Neves 2012b 


M, 


[M ] 


1.335±0.03 


1.047± 0.024 


0.7 


(Girardi 20001 


l°g£ 


[cgs] 


4.26 ±0.13 


4.28 ±0.10 


4.74 ± 0.02 


(Girardi 2000 Sousa 2011b 


age 


[Gyr] 


1.05 ± 0.64 


7.3 ± 1.2 


4.4 ± 4.0 


(Girardi 20001 


vsin i 


[km s~'] 


8.1 


2.7 


0.57 


dGlebocki 20051 CORAVEL) 


log^ K 




-4.85 ± 0.07 


-4.99 ± 0.08 






Ao,(log^ K ) 


[days] 


7 


24 




(|Noyes 1984|> 




1000. 2000. 5000. 

Period [days] 

Fig. 2. Periodogram of the RV data for HD 103774. A significant 
peak is visible at x 5.9 days. The topmost dashed line is the 
lCT 4 false-alarm probability, the lower dotted lines correspond 
to 10~ 3 , 10 -2 , and 10 _1 false-alarm probabilities. 

HD103774 haniSmean 
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Fig. 3. Keplerian fit of the radial velocity data of HD 103774 as 
a function of orbital phase. 



cant, and we see no correlation of the radial velocities with the 
activity indicators S mw (Baliunas 19951) or logiRnx) (Pearson 
coefficient < 0.1), nor with the bisectors or the FWHM. The 
time series of the activity index log(R HK ) shows strong vari- 
ations of the index with time (Fig. |4i, indicating that the star 
activity level might increase considerably. 
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Fig. 4. Variation of the radial velocity (upper panel) and the ac- 
tivity index (lower panel) of HD 103774 with time. 



Because we have enough data points, we separately analyzed 
three groups of RV data in three disjoint Julian date bins. From 
orbital fitting we have obtained the same periods for the three 
groups of data within the experimental uncertainty of w 10~ 3 
days. The same is true for the phase of the orbit, although in 
this case the uncertainties are larger (» 2 days). Therefore the 
detected period is constant within the time span of the measure- 
ment, i.e., 7.5 years. 

The period of w 5.9 days is quite close to the estimated ro- 
tation period of the star (Noy es 1984) > of » 7 days. However, the 
lack of periodicity of the bisector span, together with the lack 
of correlation between bisector and radial velocity, and the fact 
that the orbital period is strictly constant along the 7.5 years of 
data gives us confidence that the most likely cause of the radial 
velocity variation with the period of « 5.9 days is a low-mass 
companion. 

When looking at the residuals we still notice a significant 
scatter of the data. Moreover, the periodogram of the resid- 
uals shows an excess of power corresponding to a period of 
w 750 days, jointly with its one-year alias at * 250 days. We 
significantly improved the overall fit to the data by fitting two 
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Fig. 5. Activity index log(R HK ) of HD103774 as a function of 
the residuals from the single-Keplerian fit. A strong correlation 
is evident. The value of the Pearson correlation coefficient is 0.6. 



Keplerians. However, we detected a strong correlation between 
the radial velocity residuals to the single-Keplerian fit and the 
activity index log(RnK) (Fig. |5]l with the value of the Pearson 
correlation coefficient which is equal to 0.6, and found a peak in 
the periodogram of the activity index corresponding to this pe- 
riod. For these reasons we attribute the variation of the residuals 
to the one-Keplerian fit to stellar magnetic cycles dSantos 20 1 01 > . 

The highly varying log(RnK) reaches values as high as -4.7, 
lending further support to the suggestion that the activity of this 
early-type, relatively young star could be responsible for the 
large x 1 an d the scatter of the residuals around the one-planet 
fit. 



3.2. HD109271 

We measured the radial velocity of the star HD 109271 95 times 
in 7.5 years. The RV varies by 6ms~ l RMS, well in excess of the 
average RV error, which is below \ms~ l . The periodogram of 
the observed data is shown in the top panel of figure [6] and two 
peaks are clearly visible, one of which is above the 10~ 4 FAP 
level. We fit a Keplerian to this signal and obtained a good fit for 
a period of w 7.85 days. 

After fitting the first period, looking at the periodogram of 
the residuals (figure |6j middle panel) we notice a second peak, 
corresponding to the second-highest peak in the periodogram of 
the observations, above our FAP limit. This signal corresponds 
to a variation with a period of w 30.98 days. 

The phased radial velocity curves of the two planets result- 
ing from a two-Keplerians (simultaneous) fit and after subtrac- 
tion of the RV signal of the other planet are shown in figure 7] 
Neither of the signals is correlated with the bisector, with the 
activity index, or the FWHM of the CCF; we also verified that 
the periodograms of these parameters have no power at periods 
corresponding to the identified orbital periods. Moreover, both 
signals were identified using the entire data set and also using 
three disjoint subgroups of the entire data set, indicating that 
they are present during the whole time span of our observations. 
Therefore we consider the origin of the two signals to be two 
low-mass companions orbiting the star HD 109271. 

After the two-planets fit to the HD 109271 RV data, the pe- 
riodogram of the residuals shows two signals above our FAP 
threshold at Pw 198d and Pw 430d (figure [6] bottom panel). 
Within the uncertainties, these periods are the one year alias 
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Fig. 6. Periodograms of HD 109271. From top to bottom: peri- 
odogram of the data, the residuals after the one-planet fit, and 
the residuals after the two-planets fit. The dashed line marks the 
10~ 4 FAP, the dotted lines, from top to bottom, the 10 3 , 10~ 2 , 
and 10" 1 FAP. 



of each other. When constraining the analysis to consider time 
spans shorter than one year but longer than 200 days, both peaks 
disappear, and they reappear together only when the time span 
is longer than w 600 days, an indication that the signal could 
be the peak at m 430 days, and the one at m 198 days could 
be the alias. The third signal, which could be associated to a 
1.3 Neptune-mass planet on a 1 AU eccentric (e« 0.36) or- 
bit, is very noisy and at the limit of detectability by our survey 
{Kl < 2ms~ l ,< 5RV >~ lm^ 1 ). Although we cannot asso- 
ciate this signal to stellar activity, we need a larger data set and 
a higher precision to eventually confirm it as a companion to 
HD 109271. 



3.3. BD-061339 

We collected 102 radial velocity measurements of the star BD- 
061339 in about eight years . The average RV uncertainty is 
2.7 ms _I . Figure [8] shows the periodogram of the observations, 
clearly indicating an excess of power at m 125 days, well above 
the 10~ 4 FAP level. 

A Keplerian fit to the data allowed us to reconstruct an or- 
bit whose RV variation as function of the orbital phase is shown 
in Figure [9] The orbital fit is very robust, no correlation is seen 
between the radial velocities and the activity index, the bisector 
span or the FWHM of the CCF, and subsamples of the data se- 
lected in Julian date show a periodic RV variation with the same 
period within the uncertainties. 
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Fig. 7. Radial velocity curves of planets HD109271b (top panel) 
and HD 109271c (bottom panel) as a function of phase. 
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Fig. 8. Periodogram of the observations of BD-061339. The pe- 
riod at « 125 days is clearly visible well above the 1CT 4 FAP 
level. 
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Fig. 9. Orbital fit of the 125-day signal to the BD-061339 RV 
data. 



The periodogram of the residuals shows an excess at Pa: 3.9 
days with the signal just above the 10~ 4 FAP level. 

The peak in the periodogram at as 3.9 days is well defined, 
and the fit with a two-planets model converges and significantly 
improves the^ 2 . The F-test gives a 99% probability of improve- 
ment of a two-Keplerian simultaneous fit over a one-Keplerian. 
The uncertainty on the period from the fit is less than 10~ 3 
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Period [days] 
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Fig. 10. Periodogram of the residuals to the first Keplerian fit to 
the of BD-061339 RV data. A signal at Pa 3.9 days is visible 
just above the 10~ 4 FAP level. 
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Fig. 11. Evolution in time of the integrated flux in the Ca II H 
line of BD-061339. The star appears much less active during 
the last observation campaigns. The red points (darker) represent 
the amplitude of the Ca II H re-emission normalized at 6540A, 
while the green points (lighter) represent the S mw index. The 
latter one is heavily affected by the low S/N in the blue side of 
the spectrum. 



days due to the many periods covered through the entire ob- 
servations campaign. Because of the quite red color of the star 
(B — V = 1.321), the activity indicator log(RfjK) is not reliable. 
In an attempt to study the activity of this evolved star, we ob- 
served the evolution of the re-emission in the Ca II H and K lines 
with time. To quantify the re-emission we normalized the spec- 
tra to the continuum. As suggested by Mould (1976), the region 
around 6540A is only weakly affected by TiO bands. By direct 
inspection we verified that the contribution of the TiO bands at 
this wavelength for this star is indeed negligible, and that the 
width of the H a line is such that it does not affect the region we 
planned to use to estimate the continuum. The normalized inte- 
grated flux of the Ca II H line as function of time is shown in 
figure [11] It drops by a factor a 2 between JD 2454000 and JD 
2456000. In figure 1 1 we also show the time dependence of the 
S mw index: this index drops as well as a function of time, al- 
though the slope is significantly smaller. The S mw index is nor- 
malized to the continuum much closer to the re-emission line, 
but the signal-to-noise ratio (S/N) is very low in this region, 
and the substantial photon noise strongly influences the measure- 
ment. For this reason we attempted the normalization at 6540A. 

The star is becoming less active in recent times. Figure [12] 
shows the re-emission in the Ca II H line in April 2007 (JD = 
2454200) and September 2012 (JD = 2456200). The intensity of 
the line almost halved during this period. 
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Fig. 12. Ca II H re-emission from two spectra of BD-061339 
acquired in April 2007 and September 2012. The intensity of the 
line of the more recent spectrum is reduced almost by half. 
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Fig. 13. Radial velocity as function of the orbital phase of the 
short-period planet of the star BD-061339. 



Because of the variation of the activity level of the star with 
time, we divided our data set into two subsamples: one contain- 
ing data acquired before JD = 2454500 and the other being its 
complement. These subsamples contain about the same number 
of measurements, therefore our ability to identify the periodic 
variation induced by the potential second planet should be the 
same in both of them (the sampling is roughly similar). The nor- 
malized power at the peak in the periodogram goes from 3.9 
(arbitrary units) for observations taken before JD=2454500 to 
10.1 for observations taken after this date, the power increases 
for each new added data point. The period stays constant within 
one sigma in the two datasets. It is worth underlining that the 
uncertainty on the period is small also for the two separate sub- 
samples: less than 10~ 2 days. 

Adding the fact that we see no correlation between the RV 
residuals to the one-planet fit and the FWHM of the CCF or 
the bisector span, nor do we see significant power in the peri- 
odogram of these parameters around the a; 3.9 days period, we 
gain confidence that the second signal, which is becoming more 
evident as the star becomes less active, is indeed a planet. The 



RV as a function of the orbital phase is shown in figure 13 



This case is nice and instructive, and clearly shows how a 
planetary signal emerges from the noise of an active star when 
the activity level drops. It also reminds us that "active" stars are 
not necessarily always active, and that if at some time they are 
inappropriate for a planet search because of their high level of 



Fig. 14. Mass-period distribution of the RV-detected planets. The 
crosses represent all detected planets, the filled triangles the 
planets from this program, the asterisks the solar system plan- 
ets, and the filled circles the planets presented in this work. 



activity, at other times detection of planets around these stars 
may become possible. 



4. Discussion 

We have identified five planetary signals with msin(z') going 
from 0.6 to 7 times the mass of Neptune. Assuming they lie in the 
same region of the radius-mass space as the confirmed transiting 
planets (Schneider 2012 ), we can estimate their radius and there- 
fore their density. With this assumption, we obtain a maximum 
density for these newly discovered planets that is comparable to 
Neptune's or Jupiter's mean density, making them most likely 
gaseous planets. 

Observing the current detections in the mass-period diagram 



(figure 14 1, we see that the planets presented in this paper are 
close to the detection limit of our survey, of about * 3ms , and 
that two of them fall in the low-populated region of intermediate 
periods, between m 10 and m 200 days. This region is indeed 
mostly occupied by planets detected in multiple planetary sys- 
tems. It makes us wonder whether the presence of more than 
one planet breaks a certain symmetry that then allows the planet 
to occupy a range of otherwise "forbidden" orbital periods. 

The periods of the two Neptune-mass planets orbiting 
the star HD 109271 are near a 1:4 resonance, but this would 
require a thorough dynamical analysis to be confirmed. The 
planetary system orbiting the star HD 109271 is very similar to 
the HD69830 system dLovis 20061 1 : the periods of the planets 
b and c differ only by 10% and 2% respectively, while the 
planet masses differ by less than a factor of two, which is an 
indication that there seem to be some "preferred architectures" 
among planetary systems. Whether this is real or an effect of 
our limited understanding of the reality is an open question. 

The outer planet of the system orbiting the star BD-061339 
is particularly interesting because it falls right in the middle of 
the habitable zone (HZ) of its star. Following the method out- 
lined by Kane (2012), we estimate the habitable zone of the star 
to extend from 0.5 to 0.9 AU. The orbit of planet c, even consid- 
ering its relatively high eccentricity of 0.31, lies entirely within 
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Table 2. Orbital parameters for the planets. P: period, T: time at periastron passage, <?: eccentricity, a>: argument of periastron, K: 
semi-amplitude of the radial velocity curve, m sin i : projected planetary mass, a: semi-major orbital axis, V: average radial velocity, 
Mne as : number of measurements, < (5RV>: average uncertainty on the radial velocity measurements, cr (O-C) : root mean square of 
the residuals from the Keplerian fit. 



Parameter 




HD 103774b 


HD 109271b 


HD 109271c 


BD-061339b 


BD-06 1339c 


P 


[days] 


5.8881 ±0.0005 


7.8543 ± 0.0009 


30.93 ± 0.02 


3.8728 ± 0.0004 


125.94 ± 0.44 


T 


[JD-2400000] 


55675.4 ± 2.0 


55719.1 ± 3.6 


55733 ±8 


55220.50 ± 0.10 


55265.2 ± 17.6 


f 




0.09 ± 0.04 


n or i r\ ao 

0.25 ± 0.08 


0.15 ±0.09 


0.0 (fixed) 


0.31 ± 0.1 1 


OJ 


[deg] 


-42±64 


-64±20 


4± 120 


0.0 (fixed) 


41 ±55 


K 


[m s~'] 


34.3 ± 1.8 


5.6 ±0.5 


4.9 ± 0.4 


4.4 ± 0.6 


9.1 ±2.9 


m sin i 


[M Jup ] 


0.367 ± 0.022 


0.054 ± 0.004 


0.076 ± 0.007 


0.027 ± 0.004 


0.17 ±0.03 


m sin (' 


[M Earlh ] 


116 ± 7 


17 ± 1 


24 ±2 


8.5 ± 1.3 


53 ± 8 


a 


[AU] 


0.070 ± 0.001 


0.079 ± 0.001 


0.196 ±0.003 


0.0428 ± 0.0007 


0.435±0.007 


V 


[kms H ] 


-3.049±0.003 


-4.903±0.001 


-4.903± 0.001 


23.625 ± 0.003 


23.625 ± 0.003 


N 




102 


87 


87 


102 


102 


Time span 


[days] 


2734 


2683 


2683 


2955 


2955 


< 6RV> 


[ms- 1 ] 


2.6 


0.9 


0.9 


2.7 


2.7 


cr (O-C) 


[ms- 1 ] 


11.43 


2.05 


2.05 


4.3 


4.3 






15.7 ±0.6 


4.6 ± 0.3 


4.6 ± 0.3 


2.47 ± 0.23 


2.47 ± 0.23 



the HZ. Depending on the albedo (in the range 0.4 - 0.9) and on 
the atmospheric model (efficient / inefficient heat transport) the 
surface temperature of the planet can be between °C and 100 
"C. Although suggestive, this is a simplistic approach however, 
because it does not take into account that the albedo as well as 
the stellar flux vary with wavelength. Moreover, due to its mass, 
the planet is likely to be a gas planet, making conditions for life 
probably more difficult. 

5. Conclusions 

We have presented new results from our nine-year low-precision 
(2 - 3ms~ l ) survey of a volume-limited sample of the solar 
neighborhood with HARPS. We detected five exo-planets orbit- 
ing three low-activity, slowly rotating solar-type stars. Four of 
the detected planets are in two multi-planet systems. The planets 
have msin(z) ranging from 0.6 to 7 Neptune masses, the lightest 
being in the super-Earth mass range, underlining that our survey, 
albeit limited in precision because of a strategical choice, has the 
capability to detect Neptune-mass and in some cases even super- 
Earth planets. 

The star BD-06 1339 hosts two planets: in the case that 
sin(i) = 1, they would be a super-Earth, possibly in a tidally 
locked configuration due to its proximity to the host star, and 
a planet of a* three Neptune masses in the habitable zone. The 
case presented by this star is particularly instructive because it 
shows how one of the signals became more clear with time, as 
the stellar activity was diminishing. 

The star HD 103774, hosting one planet with m sin(z) roughly 
the mass of Neptune, has also shown signs of varying activity 
during the nine years of our monitoring campaign. A reminder 
that the stellar activity changes with time, and stars that at one 
time are unsuitable for high-precision RV measurements, might 
become more quiet at another time and offer the possibility to 
acquire precise RV data. 

The star HD 109271 hosts two planets with m sin(z) approxi- 
mately the mass of Neptune and with orbital periods and m sin(z') 
very similar to planets b and c around the star HD698430 
dLovis 20061 1 . This fact suggests that some planetary architec- 
tures might be "preferred" in Nature. 

Two out of the five planets we detected have periods in the 
range between 10 and 200 days, a region where planet occur- 



rence is generally low. Interestingly, most of the planets with 
periods in this range belong to multiple planetary systems. 

Owing to the small amplitude of the signals and to the prox- 
imity of their amplitudes, approximately 100 data points per sys- 
tem were necessary to achieve the current detections, which sup- 
ports ample time coverage and high-cadence monitoring as the 
key elements to access some of the most interesting planetary 
detections. 
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